Entanglement Entropy and Spectra of the One-dimensional Kugel-Khomskii Model 
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We study the quantum entanglement of the spin and orbital degrees of freedom in the one- 
dimensional Kugel-Khomskii model, which includes both gapless and gapped phases, using analytical 
techniques and exact diagonalization with up to 16 sites. We compute the entanglement entropy, 
and the entanglement spectra using a variety of partitions or "cuts" of the Hilbert space, including 
two distinct real-space cuts and a momentum-space cut. Our results show the Kugel-Khomski 
model possesses a number of new features not previously encountered in studies of the entanglement 
spectra. Notably, we find robust gaps in the entanglement spectra for both gapped and gapless 
phases with the orbital partition, and show these are not connected to each other. We observe 
the counting of the low-lying entanglement eigenvalues shows that the "virtual edge" picture which 
equates the low-energy Hamiltonian of a virtual edge, here one gapless leg of a two-leg ladder, to the 
"low-energy" entanglement Hamiltonian breaks down for this model, even though the equivalence 
has been shown to hold for similar cut in a large class of closely related models. In addition, 
we show that a momentum space cut in the gapless phase leads to qualitative differences in the 
entanglement spectrum when compared with the same cut in the gapless spin- 1/2 Heisenberg spin 
chain. We emphasize the new information content in the entanglement spectra compared to the 
entanglement entropy, and using quantum entanglement present a refined phase diagram of the 
model. Using analytical arguments, exploiting various symmetries of the model, and applying 
arguments of adiabatic continuity from two exactly solvable points of the model, we are also able 
to prove several results regarding the structure of the low-lying entanglement eigenvalues. 

PACS numbers: 75.10. Jm,75.10.Kt,03.65.Ud 
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I. INTRODUCTION 

Over the last few years, quantum information theory 
has come to play a major role in deepening our un- 
derstanding of strongly correlated quantum many-body 
systems^ The focus of these studies has largely been on 
the entanglement entropy which is a single number ob- 
tained from the reduced density matrix of a sub-system 
of the larger quantum system of interest^— However, the 
entanglement spectrum (i.e. the full set of eigenvalues of 
the reduced density matrix)^ has emerged as a powerful 
tool to study strongly correlated quantum systems due 
to its ability to reveal subtle topological effects, and its 
direct connection in some cases to the physical boundary 
excitations^^— It is remarkable that information about 
the excited states of a non-trivial quantum many-body 
system can be obtained from the reduced density ma- 
trix, which is constructed solely from the ground-state 
wavefunctioni&ii As the ground-state can be more eas- 
ily obtained than excited states via available numerical 
methods, understanding the implications of the ground- 
state quantum entanglement for the excited states is a 
central issue in quantum many-body theory, particularly 
in those systems with some type of topological order J^r— 

The entanglement spectrum has been studied ex- 
tensively in the fractional quantum Hall effect ^^r— 
spin chains ^ir—, other one-dimensional systems ; 28 ' 29 
topological superconductors ; 30 ' 31 and topological 
insulators The information contained in the en- 
tanglement spectrum depends crucially on how one 
partitions (cuts) the Hilbert space of the system into two 
portions. For example, a bond partition defined by a lo- 



cal cut across a bond of the one-dimensional Heisenberg 
spin- 1/2 chain reveals an uninteresting structure, while 
partitioning the system in momentum space reveals 
the underlying U(l) conformal field theory (CFT) and 
allows a connection to be drawn to fractional quantum 
Hall systems^ As a rule of thumb, one might expect 
local (in real-space) cuts to reveal more information in 
gapped systems, while non-local cuts (in real-space) may 
contain important information in gapless systems.— 

The one-dimensional Kugel-Khomskii model with L 
sites^ 

L 

H(X, Y) = Y,(0i ■ S i+ i + x)(n ■ f l+ i + Y), (l) 

i=l 

contains both gapped and gapless regions in its phase di- 
agram (as a function of X and Y)r2r— so we will make 
use of both local and non-local cuts (in real-space) of 
the Hilbert space to obtain a comprehensive picture of 
its entanglement properties. In Eq.([T]), Si represents 
the spin- 1/2 degree of freedom on site i, and fj is a 
two-state orbital degree of freedom on the same site. 
For general X and Y, Eq.Q possesses SU(2)xSU(2) 
symmetry, with a larger SU(4) symmetry realized for 
X = Y = 1/4. We note that Eq.([T|) admits an exact 
solution by Bethe Ansatz at the SU(4) symmetric point, 
where it is gapless ^£ and at the point X = Y = 3/4, 
where it is gapped^ 9 - and a topological "non-Haldanc" 
phase is realized^ - Phase diagrams for this model have 
been obtained via numerical methods and with field the- 
ory approaches.— ~— 

In this paper, we investigate the entanglement prop- 



2 



erties of the one-dimensional Kugcl-Khomskii model (JTJ) 
via exact diagonalization and analytical techniques. We 
begin with a study of the spin-orbital entanglement en- 
tropy using larger system sizes (L = 16) than previously 
reported (L = 8)M- For 8 sites, we find that we are able 
to reproduce the results of Ref. [5l[ , but we disagree with 
their physical interpretation. With 12 and 16 site lattices 
we obtain a better understanding of finite size effects 
present for L = 8 and find a phase diagram in overall 
agreement with other studies.— ~— 

In the remainder of the paper, we will describe the re- 
sults of our study of the entanglement spectrum of the 
one-dimensional Kugcl-Khomskii ground state with a fo- 
cus on its antiferromagnetic phases. A central goal of this 
work is to understand how a ground state wavefunction 
manifests its physical properties through the entangle- 
ment spectrum derived from a reduced density matrix. 
A density matrix is formed by tracing over degrees of 
freedom from the pure ground state density matrix. The 
degrees of freedom to be traced over are defined by a par- 
tition of the system in two parts. These partitions are 
loosely defined by cuts which represent artificial bound- 
aries between the two partitions. The different cuts that 
we have considered are the orbital or "rung" cut, the 
bond or "leg" cut and the momentum cut. The orbital 
cut is defined by tracing over the orbital degrees of free- 
dom from all the sites of the chain. The bond cut is 
defined by cutting the along two bonds of the periodic 
chain, leaving only half of the sites untraced. The mo- 
mentum cut is non-local in real space and utilizes a map- 
ping to a magnon representation of the ground state. 

One of our main results is that we find two "entan- 
glement gaps" with the rung cut. One gap seems to 
be characteristic of the gapless antiferromagnetic region 
of the phase diagram and the other separate entangle- 
ment gap is characteristic of the gapped dimerized non- 
Haldane phase. In both cases the low-lying entanglement 
spectrum is not related to the Hamiltonian of the "vir- 
tual edges"— if the system is viewed as a two-leg ladder 
(which would be a spin-1/2 Hciscnberg model). This is 
intimately related to the fact that the phases of ((T|) are 
driven by the "four-spin" interaction (Si ■ Si+i)(fi ■ fi + i) 
that pushes physics away from either the rung singlet or 
Haldane phase favored by a "rung coupling" Si ■ Ti . 50 ' 52 
Both the rung singlet and Haldane phase exhibit low- 
lying entanglement eigenvalues that mimic the physical 
excitations of a virtual edge . 26 ' 27 In addition, we have 
generalized the non-local momentum cut of Thomale et. 
ali 2 ^ to spin chains with orbital degeneracy. However, 
we do not find an entanglement gap in the entanglement 
spectrum at the gapless SU(4) point with this cut. We 
observe a new counting of the low-lying levels that has 
yet to be matched with expectations of the correspond- 
ing SU(4)i Wcss-Zumino-Wittcn conformal field theory 
as would be anticipated from the study of U(l) conformal 
field theories^ 

Our paper is organized as follows. In Section [Hi we 
present the phase diagram of the Kugcl-Khomskii model 



Eq.([T]), and identify special points of interest that we will 
focus on in our entanglement study. In Section IIII( we 
study the entanglement entropy and the entanglement 
spectrum of a partition between spin and orbital degrees 
of freedom. In Section HV1 we look at the bond or leg cut 
entanglement spectrum. In Section [Vl we consider a cut 
in momentum space, and in Section I VII we present the 
main conclusions of our paper. Several technical results 
related to the "orbital" cut and the momentum cut are 
presented in the two Appendices. 

II. PHASE DIAGRAM OF THE 
KUGEL-KHOMSKII CHAIN 

The one-dimensional Kugel-Khomskii model, Eq.([T]), 
is rich in both gapless and gapped phases with inter- 
esting symmetries 42r— In one-dimension this model can 
be thought of as cither a two-leg spin ladder (spin on 
one leg, orbitals on the other) or a one-dimensional spin- 
orbital chain at quarter filling. Physically, the Kugel- 
Khomskii model (in one, two, and three dimensions) nat- 
urally arises in the strong coupling limit (large Hubbard 
U limit) of transition metal oxides,— a class of materials 
which have recently been predicted to contain topologi- 
cal and other exotic magnetic phases under a variety of 
different conditions^ - — 
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FIG. 1. The phase diagram of the one-dimensional Kugel- 
Khomskii model given in Eq.fl]). Phase I is a fully polarized 
ferromagnetic (FM) state for both spin and orbital sectors. 
In phase II, the ground state is antiferromagnetic (AFM) for 
spin and FM for orbital degrees of freedom. In phase III, spin 
and orbital orders are reversed relative to phase II by X f-> Y . 
Phases IV and VI are gapped dimerized phases, and phase V 
is a gapless phase with non-trivial spin-orbital entanglement. 
The SU(4) point atV = F = l/4is indicated by a * and the 
exactly solvable point X — Y — 3/4 indicated by a + is in 
the non-Haldane phase with finite string order. 

The phase diagram of Eq.((T]) has been studied with 
various field theoretical and numerical methods42r— and 
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is presented in Fig. [TJ Phase I is a fully polarized ferro- 
magnetic (FM) state for both spin and orbital sectors. In 
phase II, the ground state is antiferromagnetic (AFM) for 
spin and FM for orbital degrees of freedom. In phase III, 
spin and orbital orders are reversed relative to phase II by 
X -H> Y. Phases IV and VI are dimerized AFM phases, 
and phase V is a gapless phase AFM in spin and orbital 
sectors with non-trivial spin-orbital entanglement. Thus, 
phases I, II, III, and V are gapless, while phases IV and 
VI are gapped. 

The lowest energy excitations of phase IV and VI dif- 
fer in the character of their lowest excited states which 
are magnons. Despite that, phase IV and VI are really 
the same phase, which is the non-Haldane phased At 
the boundary line separating IV and VI in FigJlJ the 
ground state is in fact fully gappe d 44 ' 47 and the transi- 
tion between IV and VI should be interpreted as a smooth 
crossover.— This unfortunately has caused some confu- 
sion in the literature 5 ^ where a distinction has been made 
between these two phases. In fact, we will show that the 
entanglement properties of the ground state are smooth 
across the combined region of IV and VI, once subtleties 
due to degeneracies and finite size effects in the dimer- 
ized phase arc taken into account. In this paper when we 
refer to phase IV, we will often mean both phase IV and 
VI, unless an explicit distinction is made. 

At X = Y = 1/4 this model contains a global 
SU(4) symmetry, is gapless, and is exactly solvable by 
Bethe ansatz4£ This SU(4) symmetry point has central 
charge c=3 and is described by and SU(4)i Wess-Zumino- 
Wittcn(WZW) theory at low-energies, which establishes 
it as a phase that is not simply the "sum" of two gap- 
less legs of a spin-orbital ladder (as the latter would have 
central charge c=2), but rather the sum of two SU(2) 2 
theories, each of which has central charge c=3/2.— In the 
two-leg ladder picture, it is the "four-spin" interaction 
(Si ■ Si+i)(fi ■ Ti+i) favoring dimerization that is responsi- 
ble for much of the interesting physics of this model . 50 ' 52 

At X = Y = 3/4 this model is gapped, but has ex- 
actly doubly degenerate ground states. As was shown by 
Kolczhuk and Mikeska,— either ground state has an ex- 
pression in the form of a matrix product state(MPS) and 
both spontaneously break the original translational sym- 
metry of the chain by exhibiting an exact dimerization. 
In the two-leg ladder picture, these ground states are the 
two inequivalent exact staggered dimer coverings of the 
ladder. The phase at this Kolczhuk and Mikcska(KM) 
point where X = Y = 3/4 is an example of a non- 
Haldane quantum spin liquid 5 ^ with a non-local order 
parameter: it has a non-zero string order paramete r 69 ' 70 
like the Haldane state,— but possess a finite dimer cor- 
relation that distinguishes it from the Haldane phase. 



III. ORBITAL OR RUNG CUT 

Turning now to entanglement physics, we begin by ex- 
ploring the properties of the reduced density matrix ob- 



tained by tracing out the orbital degrees of freedom. If 
one interprets the Hamiltonian ([1]) in the language of 
spin chains, as a two-leg ladder with spin on one leg and 
orbital on the other, the reduced density matrix would 
be obtained by tracing out one leg. Such a partition 
of degrees of freedom corresponds to a cut through the 
rungs of the ladder. To make connection with previous 
work on two leg ladder spin chains, we will refer to this 
"cut" of the Hilbert space as the "rung cut" . We will 
also sometimes refer to the orbital degrees of freedom as 
pseudospin degrees of freedom. Previously this rung cut 
was studied with exact diagonalization by Poilblancj 2 ^ 
and Lauchli and Schliemann^ 7 - for a system with only 
two-spin interactions. Their results are therefore rather 
different from ours, where the "four-spin" interaction is 
crucial to the physics. 

A. Entanglement Entropy 

In this section, we present and analyze our numeri- 
cal results based on exact diagonalization for the spin- 
orbital entanglement entropy S en t(X, Y) (obtained from 
the reduced density matrix, p rc d, with the "rung cut") 
as a function of the model parameters X and Y . The 
entanglement entropy, also known as the von-Neumann 
entropy, is obtained as 

S cnt (X,Y) = -jTr{p red {X,Y)]n[p red (X,Y)}}, (2) 

where Tr denotes the trace over the (real) spin degrees of 
freedom of the reduced density matrix and we have cho- 
sen to normalize by the chain length L. Strictly speaking 
S ent is an intensive quantity and should really be termed 
the entropy density. However, for convenience and com- 
parative purposes, we will just refer to S en t as the en- 
tropy. In all our numerical calculations we use periodic 
boundary conditions. To avoid ground state degeneracies 
in region V, we limit ourselves to chain lengths L of mul- 
tiples of 4. System sizes of L = An + 2 with n an integer 
possess additional degeneracies at the SU(4) point, that 
do not reflect the behavior in the thermodynamic limit, 
so we do not consider those system sizes 4°. 

The entanglement entropy obtained from the orbital 
or rung cut has been studied before in the pioneering 
work of Chen et alM- for L up to 12 (though only data 
and figures are presented for L = 8). Based on the sys- 
tem sizes they studied, the authors of Rcf. [5l[ concluded 
that S'ent (X, Y) can reveal phase boundaries between the 
various ground states in Fig fT] and provide a finer charac- 
terization within the ground state phases earlier found by 
Itoi et al£± A key result of Chen et al. is that S'ent (^j Y) 
reproduces the phase diagram of Zheng and Oitmaa^ 7 - 

In particular, Chen et al£^ argued (based on the en- 
tanglement entropy) that the transitions between regions 
I, II and III with V in FigQ] arc first order. The first or- 
der nature of these transitions is manifested as a large 
discontinuity in 5 en t from zero (in phases I, II, and III) 
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to a intermediate value (for region V)^ This result is 
upheld by our study, and others. The boundary between 
regions VI and IV in FigJT] also appears as a discontinu- 
ity in Sent; but to a lesser degree. However, in Ref.[5l| 
it was not appreciated that both region VI and IV are 
in the same phase and are only distinguished by the na- 
ture of their lowest lying (gapped) excitations which are 
magnonsi 46 ' 47 In fact, the "phase boundary" separating 
them has only gapped excitations and thus their ground 
states are adiabatically connected. Hence, it is surprising 
to find such a large quantitative change in 5 ent . The more 
severe transition between the gapped and gapless phase 
which corresponds to the boundary separating regions V 
and VI was identified with "ridge lines" in the topogra- 
phy of S en t(X, V)^i where a ridge is taken to mean a line 
of points where there exists at least one direction in which 
S cn t is a local maximum. It was then claimed that the 
ridges meet at the SU(4) symmetry point, X = Y = 1/4, 
which was also claimed to be a distinguished local max- 
imum. 

In this early study of the von-Neumann entropy of the 
spin-orbital entanglement in the Kugel-Khomskii chain, 
exact diagonalization data was reported almost entirely 
for chains of lengths L = 8 and to a limited degree 
L = 12j£i In our study we have reproduced their data 
for L = 8 chains and also extended our study to include 
L = 12 and L = 16 chains. Our results on larger sys- 
tem sizes show some of their conclusions are incorrect 
and their misidentified phase boundaries are due to level 
crossings between almost degenerate ground states and 
finite size effects. Degenerate ground states in phase IV 
and VI are a symptom of the spontaneously broken sym- 
metry of the non-Haldane phase caused by dimerization. 
However this degeneracy is only approximate at small 
system sizes and level crossings between these ground 
states can have significant effects on the von-Neumann 
entropy. We shall argue that with the exception of the 
phase boundaries between regions I, II, and III with V, 
Sent by itself cannot sharply reveal many of the phase 
boundaries with finite system size studies. Instead, the 
entanglement spectrum proves to be a more powerful tool 
by revealing more structure. 

For completeness, we briefly review our numerical pro- 
cedures. To compute the ground state via exact diag- 
onalization we consider the Y\ Sf = Y]- rf = sector. 
We make use of the fact that linear momentum is a con- 
served quantum number and will take the value of either 
7r or due to time reversal symmetry. In the study of 
L = 8 chains, we have performed exact diagonalizations 
in a grid of points spaced AX = AY = 0.005 apart. In 
our L = 12 study, because of computing constraints we 
were limited only to a resolution of AX = AY = 0.01. 
In our work we focus only in the region < X, Y < 1, or 
regions II, III, IV, V and VI of FigHJ 

In Figj2]Ja) we reproduce the main result of Chen et 
al£^ but with much greater resolution. The line of dis- 
continuities in S^nti which was thought to separate re- 
gions V and VI, is now made more apparent. However, 



there is also a second line of discontinuities which is re- 
vealed and which contains the exact staggered dimer KM 
point. This second line was not interpreted to correspond 
to a phase boundary previously. However, we identify 
both these lines of discontinuities as corresponding to 
level crossings between the total momentum K = and 
K = 7r ground states in the gapped phase. Field theory 
studies predict a spontaneous dimerization in the gapped 
phase which requires at the very least, doubly degener- 
ate ground states^ But with finite system sizes, this 
degeneracy is lifted and the eigenstates are only distin- 
guished by their linear momentum, in this case K = 
and K = tt. Previous numerical studies^ 4 - have revealed 
these almost degenerate lowest eigenstates in this region 
of the phase diagram. But with increasing system size 
we expect these artifacts to diminish in severity and this 
is indeed what happens with L = 12. 

Shown in Fig(2]Jb) is a surface plot of S en t(X,Y) cor- 
responding to chains of length L = 12, where S en t(X, Y) 
is now continuous as a function of X and Y with the ex- 
ception of the boundaries with the ferromagnetic phases, 
II and III (in the orbital sector and spin sector, respec- 
tively) . In place of the above mentioned two lines of dis- 
continuity for L = 8, are kinks with a discontinuous first 
derivative. Moreover, the position of these lines of kinks 
arc shifted from the that of the discontinuous jumps seen 
in the L = 8 study and there are more kinks in the L = 12 
case. This can also be clearly seen along the X = Y line 
shown in Fig 13] Thus, we see that there are significant 
qualitative and quantitative differences with the smaller 
(L = 8) system size. 

Shown in Figsj2^c-d) are contour plots of S on t(X, Y) 
with arrows drawn which represent the streamlines in 
the direction of steepest descent. These lines allow one 
to estimate the position of the ridge lines which corre- 
spond to the separatrices of the streamlines. One ridge 
line extends along the X = Y diagonal and two others 
lie almost parallel to the X and Y axes. Also highlighted 
are the SU(4) and exact staggered dimer KM point, in- 
dicated with a red and blue dot, respectively. In partic- 
ular, the SU(4) symmetry point does not seem to be a 
distinguished local maximum from the other points on 
the X = Y line in either the L = 8 or the L = 12 system. 
If one views the contours in the increasing X — Y direc- 
tion, local peaks and valleys are recognized as concave 
and convex contours respectively The point where the 
contours change their curvature is according to Ref.(5l| 
where the SU(4) point should lie. However, the SU(4) 
point as marked in Figj2^c,d) does not lie on such a point. 
Furthermore, the point at which the ridge lines meet lies 
beyond X = Y=l/4in the increasing X and Y direc- 
tion. These observations are in contradiction with the 
conclusions of Chen et al^ 

The level crossings in regions IV and VI may be iden- 
tified by a closure of the energy excitation gap between 
the two lowest eigenstates in the J2 i Sf = r? = sec- 
tor. In FigJ3]we present plots of these energy gaps along 
the X = Y line in the region of level crossings. Also 
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FIG. 2. (color online) The spin-orbit entanglement entropy S en t(X, Y). Surface plots of S an t(X, Y) for (a) L — 8 and (b) L = 12 
system sizes. The red and blue dots are indicative of the SU(4) symmetry (X = Y = 1/4) point and the Kolezhuk-Mikeska(KM) 
point (X = Y = 3/4), respectively. Note the two lines of discontinuous jumps in (a), in which the KM point (blue dot) lies 
on one of these lines. Contour plots for (c) L = 8 and (d) L = 12. The pink arrows represent streamlines in the direction of 
steepest descent. 



shown are S cn t(X,Y) curves. It is clear the positions 
of the level crossings match the discontinuous jumps or 
kinks in S ent (X,Y). However, the position of the kinks 
for the L = 12 data are harder to distinguish. Nev- 
ertheless, we can locate these points from the positions 
where the second derivative Sent is peaked and they do 
agree with the positions of the gap closures. Thus we can 
confirm that indeed the large qualitative changes in the 
entanglement entropy can be attributed to levels cross- 
ing between almost degenerate ground states. Moreover, 
each level crossing is always accompanied by a change in 
the time-reversal invariant linear momentum. 

In FigfJJa) traces of S ent (X, Y) along the X = Y line 



for all chain lengths we studied are presented and com- 
pared. With increasing chain length, the regions cor- 
responding the to gapless phase V have ever increasing 
spin-orbital entanglement entropy. Note that in this re- 
gion the L = 16 ground states are non-degenerate and 
have K = as demonstrated in Fig(4]Jb) where the 
relevant ground state energies are plotted. Conversely, 
the gapped regions IV and VI exhibit ever decreasing 
entanglement entropy with increasing system size. In 
this region, the differences between L = 16, K = and 
L = 16, K = 7r are small, both in their ground state en- 
ergy, i?o, and their entanglement entropy, Sent- This is 
again consistent with the expectation that with increas- 
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FIG. 3. (color online) (a) The excitation gap between the two 
lowest eigenstates for the L = 8 chain in the S* = r* 
U sector along the X = Y line in the region of level crossings. 
Shown for qualitative comparison is Sent in scaled units. In 
regions where the gap is non-zero, the linear momenta of the 
absolute ground state is also labeled, (b) The analogous result 
for L = 12 chains. 




FIG. 4. (color online) (a) S ea t(X,Y) along the X — Y line 
for chain lengths L = 8, 12, 16. The data presented for L 
I fi includes ground states from both time-reversal invariant 
linear momenta K = 0, n. (b) The ground state energy Eo of 
L = 16 chains from the relevant linear momenta sectors. 



There is in fact a simple relationship between the ES and 
S ent given by 



ing system size the numerical eigenstates tend to exact 
degeneracy in the gapped phase. 

Since in the thermodynamic limit of the gapped, 
dimerized non-Haldane phase, the two lowest eigenstates 
must be doubly degenerate, we do not expect the kinks 
or jumps in Sent to persist in the infinite size limit unless 
it involves a transitions between the AFM and FM phase. 
This implies a weakness in determining phase boundaries 
by solely relying on large qualitative changes in the en- 
tanglement entropy and the subtleties involved with al- 
most degenerate ground states. In the next section we 
examine the entanglement spectrum, which allows one 
to gain more insight into the nature of the ground state 
phases and the transitions between them. 



B. Orbital or Rung Cut Entanglement Spectrum 

The entanglement spcctrum(ES) is the spectrum of an 
"entanglement Hamiltonian" , H cnt , obtained from the re- 
duced density matrix via the relation^ 

p Ied = exp{-H e nt}, (3) 

and we denote the n th eigenvalue of H ent by £„. One may 
think of the ES as just the eigenvalues of the reduced 
density matrix arranged on a negative logarithmic scale. 



Sent = y ^2 £™ CX P(~£»)> ( 4 ) 
n 

where the sum runs over the entire ES. 

In the remainder of this section, we will study the ES 
derived from the orbital or rung cut and plot the spec- 
trum as a function of the momentum of the sub-system, 
which is a good quantum number due to the translational 
symmetry of the cut. In later sections, we will also study 
the entanglement spectrum for a bond cut, and a cut in 
momentum space. 

In our study of the rung cut ES, we will consider 
first the effects of symmetry on the ES and discuss the 
structure of generic spectra in the (X, Y) phase diagram. 
Then we will specialize to the case where there is addi- 
tional Z2 symmetry on the X = Y line. Next we will 
focus on a few special points on this line, such as the 
exactly solvable KM point at X = Y = 3/4 and the 
SU(4) symmetry point at X = Y = 1/4. We will study a 
few key slices of parameter space, such as along the line 
< X = Y < 1 and X + Y = 0.4 for < X < 0.4 as 
well. The latter line starts in region II of Fig[TJ and con- 
tinues through region V before ending in region III. Then, 
we will study the correlation functions derived from the 
Schmidt vectors (i.e., the eigenvectors of H cnt ) in the 
gaplcss phase V, in particular at the SU(4) point. Lastly, 
we will study the changes to the ES when an exchange 
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or rung coupling between local orbital and spin is intro- 
duced. 



1. Generic Symmetries of the Entanglement Spectrum 

We first establish the notation for our discussion on 
the ES. Let H1/2 denote the Hilbert space of a chain of 
L spin-1/2 moments. To make connection with previous 
studie d 26 ' 41 ! 49 ' 50 we will for convenience take the two leg 
ladder interpretation of the spin-orbital Kugcl-Khomskii 
chain, Eq.(TT]). From this perspective, the spin-1/2 chain 
with orbital degeneracy may be thought of as a two leg 
ladder system with one chain representing a collection 
of L spin degrees of freedom (Us) and the other chain 
representing L orbital or pseudospin degrees of freedom 
(H T )- The total combined Hilbert space is then Hst '■= 
Us ®Hr- Thus p rcd [*] = Tr w J*)(*| is the reduced 
density matrix of an orbital or rung cut for the ground 
state (GS) |#). 

Symmetries of the GS implies certain symmetries 
of the ES. It should be noted, however, that a symmetry 
of |\&) is not necessarily a symmetry of the Hamiltonian 
H(X,Y), but in the cases where \ty) is non-degenerate, 
a symmetry of H(X, Y) turns out to be a symmetry of 
|\&) modulo a phase. We shall discuss the effects of each 
symmetry on the ES. 

Firstly, translational symmetry allows Predt*] to be 
block-diagonalized by crystal momentum quantum num- 
bers if |^) itself is a momentum eigenstate. Hence the 
ES may be resolved according to momentum to reveal 
interesting structure. In fact, the ES of quantum spin- 
ladders^ 6 - and Projected Entangled Pair states (PEPs) 
on cylindrical geometries 1 ^ had previously been studied 
in this manner. 

The block diagonal structure is most easily revealed 
by expanding \^) in a Fourier basis of Hi/2(— Us,Ht)- 
The allowed momenta are p = 2-k'j- with n — —L/2 + 
1, . . . , L/2. We denote such a basis by {pa p } where {a p } 
for fixed p is an orthonormal basis for the momentum 
subspace p. With this basis we can expand any wave- 
function of %st with the tensored basis {\pa p ) <g> \q(3 q )}. 

Now let |\&) be a ground state of (JTJ with momentum 
K. Expanding in the Fourier basis gives 

I*) = E E ^W<^ bV ® kPq) ■ (5) 

p+q=K a-p.fin 
mod 27r 

Forming |\&)(\&| and then taking the trace over H r yields 
a block diagonalized reduced density matrix, 

Prodi*] 

= E { E (*P",;«S«) (%^ qPq ) \ ]pap) (pa' p \ , (6) 



with q = K — p. 



Strictly speaking, the dimerizcd ground states spon- 
taneously break translational symmetry and thus can- 
not have their ES momentum resolved. But symmet- 
ric and anti-symmetric superpositions of such degenerate 
symmetry broken states are momentum eigenstates with 
momentum K = and K = tt respectively. In a later 
section that analyzes the entanglement spectrum of the 
Kolezhuk-Mikeska(KM) point (X = Y = 3/4) these sub- 
tleties will be discussed in detail. Nevertheless, in our 
finite size exact diagonalization studies, the degeneracy 
in the AFM gapped phase IV is often inexact except per- 
haps at level crossings between K = and K = it ground 
states as was discussed in Sec lIII Al This was observed in 
Fig|2] as lines of discontinuous jumps or kinks in Sent 
which also contained the KM point(blue dot). Never- 
theless, generically the non-degenerate numerical ground 
states obtained by exact diagonalization in the gapped 
AFM region are one of the two time-reversal momentum 
groundstates with a small but finite gap between their 
energies. 

The ES can also be read off from a Schmidt decompo- 
sition (SD) of |\&) with respect to the partition between 
Us and % T . When the translational symmetry is also 
taken into account, the SD can take the following form 

i*> = EE e " e * p/2 i^>®i^-*>' ( ? ) 

where K = 0,7r due to time reversal symmetry. {4>p} 
and their duals {4> p } arc the Schmidt vectors in the mo- 
mentum sector p and (K — p) of Us and "H T respectively. 
Moreover, they each independently form an orthonormal 
basis. Here is the entanglement eigenvalue that cor- 
responds to (ftp and we denote the full ES by the set of 
pairs { (p, £0 )} which runs over all momenta and Schmidt 
vectors. For a fixed |\&) and p, a given cj) p uniquely de- 
termines (f>p. However, the above decomposition is not 
unique because of the phase degree of freedom in defining 
(j) p and <f> p . For example the transformation <fi p — > e l9 (j> p 
and <f) p — > e~ l6 <t> p leaves the SD invariant. Thus, the bar 
over <f>p should not be mistaken for complex conjugation. 

The next crucial symmetry of |\&) is spatial inver- 
sion (I) that maps sites i = 1, . . . , L according to i 1— » 
— i mod L. This is a unitary transformation and in mo- 
mentum space this maps p 1— > — p mod 2n. This symme- 
try will lead to a reflection symmetry in the momentum 
resolved ES such that {(p,^ p )} = {(—P,£<t> p )}- To sec 
this, note that for a non-degenerate GS I*), it must be 
that 2j#) = (-l)" x |*) where n x = 0,1. This follows 
from X 2 = 1 and the non-degeneracy assumption. Next 
applying X to the SD of * yields 

X \^= EE c_£ ' f/2 i(-^) ® \(k+pW p ) 

p 4> P 
= (-l)^|*). 

(8) 
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Which implies 

i*>=EE e " £ '' /J (- 1 )" I i(-j , w ® \{k+pW p )- 

v <t> P 

(9) 

Next we compare Prcdl^ 7 ] and /9 rc d[(— l) nx X\&] which are 
two alternate expressions for the reduced density matrix 

Prod [*] = E C ~^" I 

P <t> P 

p red [(-l)"-J*]= J2 E e~ U »\(-P)<l>r)((-P)U (10) 



Hence the spectrum of entanglement energy levels at p, 
{£,4> P } mu st be identical to the spectrum at — p, {£<j,_ p }- 
Thus we have a p — > — p symmetry in the entanglement 
spectrum for non-degenerate groundstates. We will de- 
note equivalences of spectra by the following suggestive 
notation, 

= {(-p,u P )} = {(PiU- P n (ii) 

where p ranges over all momenta and <p p over all the 
associated Schmidt vectors. Note that non-degeneracy 
played a crucial role in ensuring that X maps a ground 
state back to itself modulo an overall sign. 

Next, due to the SU(2)x SU(2) spin and pseudospin 
(orbital) rotational symmetry, a GS \^f) if non-degenerate 
must necessarily be a singlet in both spin and pseudospin. 
That is Sl ot = r t 2 ot = 0. Tracing over T-L T keeps p ro d in 
the singlet Sf ot = sector of rls- Thus, the momen- 
tum resolved entanglement eigenvalues of p re d are ex- 
pected to be generically non-degenerate unless additional 
symmetries or accidental degeneracies dictate otherwise. 
That is, there are no degeneracies to be associated with 
spin and pseudospin rotations. This should be contrasted 
with the ES derived from the gapped two-leg Hciscnberg 
ladder^ There the GS instead satisfies the less restric- 
tive (Stot + ^tot) 2 = condition. This then leads to the 
ES having support in more than one S^ ot sector of TLs- 
Here, the orbital r-chain is identified with the second leg 
of the ladder which is traced over. 

Lastly, the S^ ot = condition dictates the maximum 
total number of physically significant entanglement en- 
ergy levels from our numerically computed spectra. Let 



-^singlet denote the dimension of the S^ ot = subspace of 
rls which grows with L. /V s i ng i e t may be calculated by 
counting the number of standard Young tableaux with 
shape [L/2,L/2] or 2 rows and L/2 columns. The ES re- 
veals its physically relevant structure with only the low- 
est iV s i n giet levels. In practice our computations generate 
more eigenvalues than Af s i ng i et because we work in the 
more convenient S* ot = r t z ot = basis. These additional 
eigenvalues of p IC( \ should be zero in value based on phys- 
ical grounds. But in our numerics they are seen to be 
finite but exponentially small compared the the rest of 
the spectrum. 
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(a) X=0.22,Y=0.17 
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(b) X=0.76,Y=0.63 
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FIG. 5. (color online) Orbital or rung cut entanglement spec- 
tra against linear momentum p at generic points in the AFM 
phases for L = 12. The o symbols denote singlet levels, 
(a) (X,Y) = (0.22,0.17) from the gapless V phase. There 
are 32 low energy states across all p that are separated by 
the gap A£hi g h from a continuum of high energy states, (b) 
(X, Y) = (0.76, 0.63) from the gapped VI phase. There are 
only 2 low energy levels at p — 0, n which are separated by 
a gap A£i ow from a continuum of high energy states. The 
position of A£hi g h is above that of A£i ow . 



2. Generic Entanglement Spectra 

Our calculations of the orbital or rung cut ES from 
exact diagonalization of generic points in the AFM re- 
gion do indeed exhibit the symmetry properties of the 
last subsection. Shown in Figj5] are entanglement spec- 
tra taken from a L — 12 chain derived from ground states 
deep in the gapless and gapped AFM phases. Both spec- 
tra are p — > —p symmetric and are entirely composed of 
non-degenerate singlets. However the ES is interesting 
in other ways. The striking feature is an entanglement 
gap across all momenta exhibited by both spectra. In the 
gapless case FigJSJa), this gap labeled by A£hi g h separates 
a set of 32 low lying energy levels from a continuum of 
high energy levels. While in Figj5^b) a different entangle- 
ment gap A£i ow separates a set of high levels from only 
two distinguished lowest levels at p = 0,tt. The corre- 
sponding ES for L = 8 (L = 16) is qualitatively similar 
but with less (more) entanglement levels. In phase V, 
data from L = 8,12 and 16 reveal an empirical relation 
2-L/2-1 j? or nurn ^) er f lowest lying levels which are 
distinguished by the entanglement gap A£high- In phase 
IV, there are still only two lowest levels at p = 0, tt. 

Remarkably the rung cut ES of FigJS] does not in 
any way resemble the gapped two-leg Heisenberg spin- 
ladder^ even in the gapped IV phase. The ES does not 
seem to exhibit any similarities, quantitative or qualita- 
tive with the real spectrum of a single Heisenberg chain. 
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More specifically, triplets and higher degenerate multi- 
ples are absent from the spectrum. This could be in- 
terpreted in two different ways. On the one, hand it 
could mean the failure of the conjecture that relates the 
entanglement spectrum to the real energy spectrum of 
the (virtual) edges.— On the the other hand, it could 
be that the conjecture still holds true but the boundary 
or edge Hamiltonian is not the single AFM Hciscnberg 
chain as one would expect by decoupling the spin and 
orbital chains in (Q]). 

The lack of triplets and higher degenerate multiplets 
can be traced back to the restriction to the S 2 ot = sec- 
tor of the entanglement eigenvectors of p re( j as was shown 
in the previous section with symmetry based arguments. 
Hence we believe the second scenario to be more likely 
and the special nature of the Hamiltonian with the four 
spin interaction is the cause for the complicated bound- 
ary Hamiltonian. 
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(a) X=0.05,Y=0.25 
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(b) X=0.25,Y=0.05 
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3. 1^2 Symmetry of the Entanglement Spectrum 

In this subsection we consider the effects of the Z2 
symmetry operator which exchanges spin (S) and orbital 
(r) degrees of freedom. We denote such a linear operator 
by T which acts on Ust by 



20 r 



(c) X=0.62,Y=0.62 



(12) 



for all a e Hs and (3 £ H T - Moreover, T is unitary 
and F 2 = 1. Effectively, J- exchanges X <h> Y when act- 
ing upon H{X,Y). That is TH{X,Y)T = H(Y,X). 
Thus a ground state |\&) of H(X,Y) is mapped to a 
ground state of H(Y, X) under T . Applying T to the 
Schmidt decomposition of |\E') in ([7]) yields the symme- 
try S en t(X, Y) = S cn t(Y, X) which was noted by Ref. 51 1. 
However, it is not necessarily true that the ES derived 
from the GS of H(X, Y) is identical to the one derived 
from H(Y, X). It is only that true that both spectra pro- 
duce the same von- Neumann entropy Sent- In fact, using 
symmetry arguments similar to that of Section IIIIB II it 
can be shown that when K = ir they are shifted by tt 
relative to one and another, otherwise they are identical 
when K = 0. These arguments arc presented in Ap- 
pendix [3] Shown in Fig|6ja-b) are entanglement spectra 
exhibiting this relative shift. 

Next, we specialize to the Z2 symmetric X = Y line 
where J 7 is a symmetry of H(X,X). In this instance 
there may be double degenerate levels (doublets) present 
in the momentum resolved ES. These doublets occur at 
p = ±7r/2 but only when K = n. Shown in Figjf^c- 
d) arc entanglement spectra taken from phase IV on the 
X = Y line. Fig|6fc) comes from a GS with K = whilst 
Fig|6|d) from one with K = it. Details of this argu- 
ment for this degeneracy are also presented in Appendix 
IA1 where it also requires that satisfy J r | v I / ) = — |^). 
These results emphasize the importance the ground state 
momentum K may have on the entanglement spectrum 
at finite systems sizes. It is also one more example of 
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(d) X=0.58,Y=0.58 
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FIG. 6. (color online) Momentum resolved entanglement 
spectra for L = 12 chains at points related by Z2 spin- 
orbital symmetry. The o symbols denote singlet levels while 
the e symbols denote doublets, (a) (X, Y) = (0.05,0.25) 
in phase V. (b) (X, Y) = (0.25, 0.05) in phase V which is 
shifted by p -> p + tt relative to (a), (c) (X, Y) = (0.62, 0.62) 
in phase IV with ground state momentum K = 0. (d) 
(X, Y) = (0.58, 0.58) in phase IV with K = tt. Note the 
7r shift and doublets at p — ±7r/2 are present regardless of 
whether or not the ground state is in phase V or IV. 



how the entanglement spectrum is much more sensitive 
to the details of the ground state than the von-Neumann 
entropy. 

In the next sections we will present our results and 
analyses on the entanglement spectra at the two special 
points along the X = Y symmetric line. These are the 
SU(4) symmetric point (X = Y = 1/4) and the exactly 
solvable Kolezhuk-Mikeska(KM) point (V = Y = 3/4). 
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4- Entanglement Spectrum for the Exactly Solvable 
Kolezhuk-Mikeska point X — Y — 3/4 

Recall that at X = Y = 3/4 an exact ground state 
of H(X, Y) in terms of a matrix product state (MPS) 
was found by Kolezhuk-Mikcska4£ The ground state at 
the KM point is exactly doubly degenerate even for finite 
system sizes. These two MPS ground states are the two 
possible staggered dimer coverings of the system, with 
no dimers between the S and r chains. Moreover, they 
each individually break the translational symmetry of the 
system. Thus, their ES may not be momentum resolved 
as was discussed in Section IIIIB II In addition, these 
wavefunctions are exact product states of simpler spin 
and orbital wavefunctions. Hence they are expected to 
yield no entanglement whatsoever under a rung or orbital 
trace. These ground states are also examples of a non- 
Haldane gapped spin-liquid with string order , i 50 

Nevertheless, by taking symmetric and antisymmetric 
superpositions of these exactly staggered dimers, momen- 
tum eigenstates with K = and K = tt which have the 
full symmetries of if (3/4, 3/4) are recovered. In fact, the 
numerical ground state wavefunctions computed around 
the KM point are always momentum eigenstates and are 
thus adiabatically connected to one of these momentum 
eigenstates. Interestingly, the entanglement is non-trivial 
for these K = 0, tt eigenstates at the KM point. Shown 
in Figf7| is the ES at the KM point for K = tt which 
is extraordinarily simple because it consists of only two 
levels with £ = In 2 at momenta p = 0, tt regardless of 
system size L. When K = there are still only two lev- 
els at p = 0,7r but these levels are shifted slightly from 
In 2 by an L dependent amount which vanishes in the 
thermodynamic limit. It thus appears that at the KM 
point there is an infinite entanglement gap (A£i ow = oo) 
which becomes finite when moving away from this point 
as was observed in generic spectra in phase IV. 

The entanglement levels at the KM point for K = 
0, tt can be computed exactly using the staggered dimer 
wavefunctions. This calculation is outlined in the rest 
of this subsection. We denote a dimer between sites i 
and j on a given chain by Dij with the following sign 
convention 



1 



|Aj> -j= (| T 4i> - \UW ■ 



(13) 



For sites 1,2,..., 2N = L there is canonical dimer cover- 
ing of the chain which we denote by 



DW^> := \D 12 D 3 4...D 2 n-i,2n) 
and its translated counterpart 



l,2i 



(14) 



D ( 



2)\ ._ 



= D 23 D 4 



45 • ■ • D 2 Na) — 



N-l 



J D2i,2i+lD 2 NS 



(15) 
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FIG. 7. (color online) Entanglement spectrum of an L — 
12 system from the exactly solvable Kolezhuk-Mikeska(KM) 
point X — Y — 3/4 for K = n. This point is in the non- 
Haldane phase with a finite string order parameter.— There 
are only two levels £ = In 2 at p — 0, tt. The K = spectrum 
(not shown) is almost identical except the levels are shifted 
slightly from In 2. Away from the KM point, the spectrum be- 
comes "dressed" with high energy levels as shown in Fig[5jb) 
where there is a finite entanglement gap A£i ow that separates 
the two sets of energies. 



The overlap between the two distinct dimer coverings is 
non-zero and can be shown to be 



D (i)| D (2)\ 2 (-l/2) 



(16) 



Pictorially, this says when the two coverings are over- 
laid on top of one another to form an unbroken chain of 
length 2N, then the non-zero overlap is 2(— l/2) N , with 
our convention. This is a special case of a more general 
formula for general dimer coverings of a lattice. Note also 
that as N — > oo, the overlap vanishes but is never exactly 
zero for finite N. By taking symmetric and antisymmet- 
ric superpositions and normalizing, we can form K = 0, tt 
momentum eigenstates from the two inequivalent exact 
staggered dimer wavefunctions: 



l*o> 



2 N-1 

^2(4^-1 + 1) 

x (|D«)® |D^) + |D( 2 >)(g> |D«)) , 

2 N-l 



n/2(4^-! - 1) 
x (|D«) ® |D^) - |D< 2 >) (8 ID^ 1 



(17) 



Their respective reduced density matrices can then be 
computed to arrive at an N dependent result for p tc d [^o] 
but a simple expression for Predt^ir] 1 ^ 
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Prcd[* ] = 

1 / 3 + 4 ^-! 2(-l) N V4, N - 1 - l\ 

2(4^"! + 1) \2(-l) N V^^ T 4 JV " 1 -1 J' 

Acd[*.] = x / 2 ) • (18) 



L = 2N 


Numerical results 


Computed from Prod^o] 


8 


0.473085, 0.975714 


0.473085, 0.975714 


12 


0.63258, 0.757621 


0.63258, 0.757621 


16 


0.677645, 0.708896 


0.677644, 0.708895 



TABLE I. Comparison between analytical result for £ and the 
numerics for the K — groundstate at X = Y = 3/4. 



Thus, for a K = ir ground state at X = Y = 3/4, 
we expect that 1/2 = will produce two degenerate 
entanglement energies with value In 2. This is indeed ob- 
served in our numerics for L = 2N = 8, 12, 16. With the 
above expression for jO r cd ['J'o] j the entanglement energies 
can be determined and compared with our exact diago- 
nalization results. These values are shown in Tabic U and 
are in almost perfect agreement. 

Lastly it is clear that as L — > oo, the overlap (fTo| tends 
to zero and the differences in the ES between \^q) and 
I^tt) should become negligible. 

5. Entanglement Spectra for the SU(4) point 

In this subsection we focus on the rung cut ES at the 
SU(4) symmetric point at X = Y = 1/4. Shown in 
FigiHa-c) are spectra taken from ground states at that 
point for sizes L = 8, 12 and 16. A noteworthy feature 
shared by the spectra is the entanglement gap A£high 
around 6 < £ < 10 across all momenta p, and the man- 
ner in which the gap becomes more well defined with in- 
creasing system size. Moreover a simple finite size scaling 
analysis of the entanglement gap taken at p = plotted 
in FigJBJd) suggest that the gap remains open in the ther- 
modynamic limit. 

Several features are worth pointing out regarding the 
systematic changes in the ES with increasing L. Firstly, 
the L = 8 ES of FiglUJa) is severely limited by the finite 
system size effects. Recall that -/V s i ng i t, the dimension of 
the St ot = subspace of "Hs, sets the maximum num- 
ber of physically relevant levels of the ES. For L = 8, 
-^singlet = 14 and this is the total number of levels ob- 
served in Fig(5Ja). Of these 14, 8 levels make up the 
set of entanglement levels below the gap. This sparsity 
of the non-zero eigenvalues of p re( j highlights the limita- 
tions in attempting to draw strong conclusions about the 
thermodynamic limit using only L = 8 data. 



Secondly, upon increasing L to 12 and 16 the spec- 
trum starts to fill up with levels above and below the 
gap. At L — 12, doublets are observed at p = ±n 
which can be attributed to the properties of a K = ir 
ground state as discussed in Appendix More impor- 
tantly, however, this filling of levels is observed to be 
systematic. In all sizes considered, the total number of 
levels below the entanglement gap follows the empirical 
relation 2 i / 2 ~ 1 . In a later section we will discuss the 
physical significance of the A£high entanglement gap for 
correlations functions and present some conjectures as to 
the origin of the 2 L / 2_1 relation. 

Lastly, as was already revealed in Fig(5ja) and 
FigJS{a,b), the A£hi g h entanglement gap seems to be a 
generic property of ground states deep in the gapless 
AFM phase V, where the 2 L / 2_1 counting is also ob- 
served at these points. Viewed in this light, the SU(4) 
point ES appears to be a rather generic ES of the phase 
V. This is rather surprising given that this point is 
known to be a critical point of a Kosterlitz-Thouless(KT) 
transition ! 43 ! 45 In fact, at points X = Y > 1/4 beyond 
the SU(4) point, the ES still yields a robust A£high en- 
tanglement gap. This fact unfortunately has not allowed 
us to conclusively identify the SU(4) point as being a 
critical point of the KT transition using only our present 
amount of ES data. Nevertheless, we believe that the en- 
tanglement gap is a distinguishing property of the phase 

V and our interpretation of our finite system size data 
requires taking into account the complexities of the KT 
transition. In the next few sections we will present our 
data and analyses on the evolution of the ES along sev- 
eral lines in the phase diagram with a focus on the region 

V to IV phase transition. 



6. Entanglement Spectra Along the X = Y Line 

In this subsection we investigate the evolution of the 
entanglement spectrum along the Z2 symmetric line < 
X = Y < 1 for sizes L = 8, 12 and 16. Shown in FigHare 
the non-momentum resolved entanglement spectra along 
this line and overlaid with 5 en t data. The data for L = 8 
and L = 12 is presented in Figj9^a) and FigOHb). For 
these plots, only absolute ground states are used to com- 
pute the ES. The spectra are also labeled according to 
their momentum K = 0, 71". Our L = 16 data set is unfor- 
tunately more limited due to computing limitations. In 
this case we have decided instead to show separately the 
spectra taken from ground states of definite momentum 
K sectors. These results arc shown in FigOUc) (K = 0) 
and FigEJd) (K = tt). As is shown in Figg^b), the 
K = ground state is often the absolute ground state 
for most of the < X = Y < 1 line. 

We remark first on the L = 8 and L = 12 results of 
Fig Ma, b). In Figj9lja) where L = 8 the total number 
of levels is limited to N s i ng \ ct = 14 whereas in Figl^b) 
where L = 12 the spectrum is more densely populated. 
Both plots very clearly reveal the level crossings between 
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(a) L=8 (b) L=12 (c) L=16 (d)Gapsatp=0 
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FIG. 8. (color online) Entanglement spectra at the X = Y = 1/4 SU(4) point at various system sizes. The o symbols denote 
singlet levels while the e symbols denote doublets, (a) L — 8. (b) L = 12. (c) L = 16. The doublets at p = ±7r/2 for L — 12 
are due to the ground state momentum K = 71". In all cases the total number of low lying energies below the entanglement 
gap obey the empirical relation 2 Z,//2_1 . (d) A finite size scaling analysis of the entanglement gap A£hi g h taken at p = 0. The 
positive intercept suggests that the gap remains open in the thermodynamic limit. 



K = and K = -k ground states. These are easily identi- 
fied with discontinuous changes in the spectrum and are 
in exact correspondence with the positions of the non- 
smooth changes in the entanglement entropy S cn t- Oth- 
erwise, away from the level crossings, the evolution of 
the entanglement spectrum is continuous. This is again 
another demonstration of the advantage of studying the 
ES over the entropy when trying to detect changes in the 
ground state. In Figsj9]^a,b) the high entanglement gap 
A£high is well defined and continuously evolves far beyond 
the SU(4) point X = Y = 1/4 until a level crossing is en- 
countered. The low entanglement gap A£i ow by contrast 
starts to develop around X = Y ~ 0.2. Eventually it 
diverges at the exactly solvable KM point X = Y = 3/4 
where the ES is composed of just two levels. In addition, 
there are also parameter regions where both A£high and 
A£i ow are simultaneously non-zero. Lastly, we note that 
generally a large A£high corresponds to a large Sent- 

Next, data on the ES for L = 16 shown in Figs|9{c,d) is 
significantly more sparse. Nevertheless in FigJS^c) for the 
case where K = 0, we have have computed many more 
data points in the intermediate region between phase V 
and IV with a focus on the phase transition. The ES 
from this set of ground states is continuous over the en- 
tire range of X = Y. It is observed that the two low- 
est levels at p = 0, ir from the KM point continuously 
evolve into the two lowest levels of the 2 L / 2 ~ 1 = 128 set 
of energies below A£high in the phase V. In FigOUc) the 
gap A£hi g h is observed to close around X = Y ~ 0.5 
which is before the KM point but beyond the SU(4) 
point. At this point, the entropy S ent is also seen to 
decrease rapidly with increasing X = Y. Based on this 
gap closure which occurs for L = 16 but not smaller sizes, 



we conjecture that with increasing L, the closure point 
of A£hi g h will tend towards smaller X = Y values. It 
is known that the real energy excitation gap closes only 
exponentially slowly when approaching the KT critical 
point from phase IV42r— This makes it difficult to pin- 
point the critical point of the transition without having 
to resort to a renormalization group analysis. Our nu- 
merics reflect this reality too, as determining the critical 
point by the opening or closure of cither A£i ow and/or 
A£high does not lead to a satisfactory determination of 
the phase boundary. Perhaps future work utilizing Den- 
sity Matrix Renormalization Group (DMRG) to access 
larger L sizes with appropriate finite size scaling analy- 
ses respecting the peculiarities of the KT transition will 
improve on this. 

Finally, we should mention that A£high is also well de- 
fined for the K = ir ground state (Figj9jd)) in the range 
< X — Y < 1/4. However these cigenstates are excited 
states for those parameters and interestingly, the set of 
low energy levels it defines actually exceeds the 2 L / 2_1 
limit. 

Recent work has numerically substantiated the claim 
that one can adiabatically connect two phases as long as 
an entanglement gap remains openj^i If that conjecture 
is true, one would conclude that the A£i ow gap being 
lower in position in the spectrum and thus more signif- 
icant than A£high, should open there and remain open 
throughout regions IV and VI. Again due to difficulties 
of the infinite order KT transition and the local spin- 
orbital degrees of freedom per site, this conjecture is 
hard to verify via exact diagonalization for the system 
sizes we could handle. Nevertheless, the continuous and 
open A£i ow gap (for fixed K) throughout region IV and 
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(c) L = 16, K = (d) L = 16, K = it 
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FIG. 9. (color online) The (non-momentum resolved) entanglement spectrum along the X — Y symmetric line for sizes L = 8, 12 
and 16. The blue (dark) points correspond to spectra from ground states with momentum K = and the green (light) to 
K = n. In (a) and (b) the absolute numerical ground state is used to derive an ES. Due to computing limitations, our data 
points for L = 16 are more sparse and we have plotted separately the ES for K = and K = ir sector ground states in (c) 
and (d). In phase V the absolute ground state is mostly the K — one (See Fig|4jb)). Overlaid as the red (solid) line is 
S en t(X, X) in arbitrary units. Marked out in (a) is the high entanglement gap A£hi g h that separates the set of 2 i//2_1 low 
energy levels from the rest of the entanglement spectrum. Also marked in (a) is the low entanglement gap A£i ow that separates 
the 2 distinguished lowest levels of the gapped AFM phase from the rest of the spectrum. The region where the entanglement 
entropy rapidly drops lies in the region where the transition from the gapless region V into the gapped regions IV and VI in 
Fig[T] occurs. 



VI provides further evidence that the phase diagram in 
Ref.fo"l| contained an unphysical phase boundary. Based 
on this we can again conclude with our ES results that 
the ground states of IV and VI are adiabatically con- 
nected once degeneracy due to dimerization is removed 
by fixing K. 



7. Entanglement Spectra Along the Line X + Y = 0.4 

We next consider the behavior of the (non-momentum 
resolved) entanglement spectra along the line X + Y = 
4/10 for < X < 0.4 which is perpendicular to the 
symmetric X — Y line. This line begins in region II of 



Figffland slices through through region V before finishing 
in region III. Shown in Fig llOl are the non-momentum 
resolved entanglement spectra along the line X + Y = 
0.4. The spectra plotted are taken from the ground state 
for the system size L = 16. Results for other sizes are 
qualitatively similar. 

For X < 0.06 and X > 0.36 the structure of the en- 
tanglement spectra collapses to a singular level at £ = 
with Sent = because these states have trivial entan- 
glement. (They are in region II and III, respectively, of 
FigfU) While for 0.06 < X < 0.36, entanglement gaps 
A£high and to a lesser extent A£i ow are seen. This path 
corresponds to a slice through region V. A£hi g h is largest 
in the "middle" of region V and appears to close at the 
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Spectra of Y+X=0.4 
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FIG. 10. (color online) The (non-momentum resolved) entan- 
glement spectrum along the X + Y = 0.4 line for L = 16. The 
red (solid) line is the entanglement entropy. For X < 0.06 and 
X > 0.36 the structure of the entanglement spectra collapses 
as these states have trivial entanglement. (They are the FM 
phases II and III of Fig.l.) For 0.06 < X < 0.36 the entangle- 
ment gaps A£high and A£i ow are seen. This path corresponds 
to a slice through region V. 



phase boundaries. Thus, as claimed in Scc lIII B"2l it ap- 
pears that the A£high gap in the entanglement spectra 
with the rung cut persists throughout the gaplcss region 
V. (Recall that the SU(4) point sits on the boundary of 
region V in the thermodynamic limit.) 



8. Correlation Functions and the A£high Entanglement Gap 

In this subsection, wc explore the physical significance 
of the entanglement gap A£high for ground state(GS) cor- 
relation functions. The greatest consequence of a large 
well defined A£high is the ability for certain GS correla- 
tion functions to be well approximated by using only the 
N(L) := 2 i / 2 ~ 1 Schmidt vectors below A£ high . Recall 
that tracing over rl T to form p re( j does not remove infor- 
mation from observables that only depend on the spin(5) 
degrees of freedom. In fact, it is these correlations that 
can be well approximated. 

To be more specific, let 0(xi, . . . , x m ) be an operator 
depending only on S X1 , . . . , S Xm . Then its GS expecta- 
tion value can be computed from 

(*|6(xi, . . . , x m )\V) = Tr Hs {0(x-L,. . . , x m )p Ted [^]}. 

(19) 

Let the set of entanglement levels {£ ra } with their cor- 
responding set of Schmidt vectors {4> n } be organized in 
ascending order £i < £ 2 < ■ • ■ < £/v 3ingl0 f Using this 
cigenbasis of p ro d, the GS expectation can be computed 



with 

^singlet 

(*|d(xi,...,a; m )|*> = e~ U (4>n\6{x u ...,x m )\(j) n ). 

71 = 1 

(20) 

Now, a well defined gap A£high separates the ES into low 
and high subsets. Moreover, the levels above A£high will 
have Schmidt weights that are exponentially smaller than 
those below it. This can be exploited to approximate (O) 
by, 

2 L/2-l 

mO(x 1 ,...,X m )\^} « ^ e-^{4> n \6{x 1 ,...,X m )\4> n ) 

71=1 

(21) 

with exponentially small corrections of order 
0( e -(£w(i)+ A W)). That is, with a large A£ high , 
one only requires N(L) = 2 L / 2 ~ X Schmidt vectors to 
reproduce all the GS correlation functions involving only 
spin(S). Although N(L) still grows exponentially with 
L, it is still slower than iV s i ng i et . Moreover, these GS 
correlation functions should obey the predictions of the 
SU(2)2 WZW field theory that describes the low energy 
dynamics of the spin sectorj 43 i 45 i 46 Phrased another 
way, the Gibbs ensemble^ /Oredf'J'] of size N(L) with the 
degrees of freedom of a spin-1/2 chain of length L can 
simulate well all the GS correlators of a SU(2) 2 WZW 
CFT, provided A£hi g h is large enough. 

Next, we discuss our numerical results on the GS cor- 
relation functions calculated from the Schmidt vectors 
or equivalently the eigenvectors of /9 re d[^]- Shown in 
Fig lllf b.c) are spin-spin (S± ■ Si+i) and dimer-dimcr 
(Si ■ S , 2)(S , j+2 • Sj+3) ground state expectation values 
for a L = 16 chain at the SU(4) symmetric point. Al- 
though we only present data for the SU(4) point, our 
results remain qualitatively similar at other points deep 
in phase V. We focus first on the open circles o which 
are correlation functions computed using Eq.pTjl. Wc 
have checked that these correlation functions are indis- 
tinguishable from those obtained using Eq. (|20|) where all 
the Schmidt vectors are utilized. Thus, we have con- 
firmed the validity of the approximation of Eq. (|2ip . 

In FigfTTTb) a decay is observed in the spin-spin cor- 
relation function which is unfortunately limited by (pe- 
riodic) boundary condition effects. We also observed pe- 
riod 4 oscillations which corresponds to a 7r/2 peak in the 
structure factor. These oscillations are consistent with 
previous DMRG calculations 4 ^ and field theory studies4£ 
However, due to the limited size L = 16, our data does 
not allow for a fit to the 3/2 power law decay predicted by 
the SU(2) 2 WZW field theory, 4 ^ Nevertheless, the pro- 
file is in qualitative agreement with Ref. [40j | and gives us 
confidence that we can expect to see the algebraic decay 
with longer chains using only the approximation (|21[) . In 
Fig lllf c) the dimer-dimer correlation function does not 
seem to display a strong decay like in the spin-spin case. 
Rather, it almost appears to be long-range ordered. But 
we believe this to be a finite-size effect because dimer- 
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(a) X=Y=0.25 



(b) Spin -spin Correlation 




FIG. 11. (color online) Correlations functions taken from the 
SU(4) symmetric point ground state for L — 16. (a) The 
entanglement spectrum where the blue o symbols denote the 
iV(16) = 128 levels below the entanglement gap A£hi g h- The 
two lowest levels with the greatest Schmidt weight are marked 
by the solid red • symbols. These levels are adiabatically 
connected to the two levels of the KM spectrum, (b) The 
spin-spin ground state correlation function, (c) The dimer- 
dimer ground state correlation function. In (b) and (c) the 
open circles o are computed using only levels below A£ h i g h 
while the solid • use only the two lowest levels of (a). 



ization is not expected to occur at the SU(4) symmetry 
point. 

We also considered the contributions to the GS cor- 
relators from just the two lowest entanglement energy 
levels. That is, we calculate (|20[) with just the first 
two Schmidt vectors (f>i and fa which have the great- 
est Schmidt weight. The results of these computations as 
shown as the solid • circles in Fig lllf b.c). These two lev- 
els are special in that they are smoothly connected to the 
entanglement levels of the KM point as shown in FigJSfc). 
Most interestingly, the resulting correlator in Fig lTTT c) 
demonstrates that their contribution dimer-dimer corre- 
lation function is surprisingly uniform and strongly long- 
range ordered. We have also computed the same correla- 
tion function at the KM point and a qualitatively similar 
profile is obtained except there the negative values at 
even separations is heavily suppressed. Also, the spin- 
spin correlator of the KM point is ultra short-ranged, 
decaying to zero beyond nearest neighbors. Based on 
these observations, we interpret the role of the remaining 
iV(16) — 2 = 126 Schmidt vectors with higher entangle- 
ment energies, as providing the necessary critical alge- 
braic spin-spin fluctuations that destroy this long-range 
dimerization order. This interpretation is also based on 
the expectation that the SU(2)2 WZW low energy effec- 
tive description describes algebraic spin-liquid behavior 
without dimerization. This also means that A£i ow is a 



better indicator of the gapped non-Haldane AFM phase 
IV, since it separates the two lowest KM levels and the 
next N(L) — 2 levels. It is also worth repeating here 
that in certain excited states in phase V, the A£ n igh can 
sometimes be observed, but in those circumstances the 
number of levels below it exceeds N(L). This suggests 
that with excited states which have more complicated 
correlations, there is an increased amount of complex- 
ity which manifests as an increased number of Schmidt 
vectors below A£high- 

We still lack any explanation, even if only intuitive of 
the A£high entanglement gap and the N(L) = 2 L / 2_1 lev- 
els below it. This is in contrast to the A£i ow gap which 
can been understood in terms of pcrturbative corrections 
to the KM ground state which makes the otherwise in- 
finite gap, finite. It would be very interesting to find a 
ground state which is adiabatically connected to region V 
where A£high = oo such that Eq. (|2"Tj) becomes exact. We 
can expect such a ground state to describe the SU(2) 2 
WZW model ground state more faithfully. Curiously, 
the number N(L) may have a combinatorical aspect to 
it. We noticed the following combinatorical re-expression 
of N(L), 



N(L) = 




(22) 



From this point of view, N(L) is then the number of ways 
to select an even number of L/2 objects. L/2 is also the 
number of nearest neighbor links or possible dimers. We 
conjecture that this relation may be of some significance 
towards explaining the A£high gap, although this is highly 
speculative at the moment. 

In most systems studied to date, an entanglement gap 
(with a different cut from the one considered here) has 
been shown to separate low-lying states with universal 
topological properties from higher states that contain in- 
formation about excitations ! 10 ' 12 " — In these cases, the 
systems have a gapped bulk excitation spectrum. It is 
thus somewhat unusual to find an entanglement gap in 
a gapless system. However, the one-dimensional spin- 
1/2 Heisenberg model with a cut in momentum space is 
an important example of a gapless system with an en- 
tanglement gap.— In that case^ it was shown that the 
gap could be related to a fractional quantum Hall state 
of the Laughlin type, as they both share an underlying 
U(l) conformal field theory. 

In Sec|Vj we study the entanglement spectra of the 
model ([I]) at the SU(4) point with a momentum space 
cut. Quite interestingly we do not find a gap with this 
cut, but we are limited by computational resources to 
only L = 8,12 and there are many subtleties associated 
with finite-size effects, as we noted in Sec lIIII A. Never- 
theless, it would be interesting to see if the entanglement 
spectrum from the rung cut and momentum cut could 
still be related to a quantum Hall system derived from 
an underlying SU(4)i conformal field theory, in an anal- 
ogous way to the U(l) conformal theory of the spin-1/2 
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Heisenbcrg spin chain connecting it to the Laughlin quan- 
tum Hall states.— 



(a) X=Y=0.75,/ rung =0.001 (b) X= Y=0.75,/ rung =0.01 

20 1 2 a 1 20 1 » — >» ^ 
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9. Perturbing the Kugel-Khomskii Hamiltonian {TJ with a 
Rung Coupling 

In a number of states with gapped bulk excitations and 
some type of topologically protected boundary excita- 
tions, such as fractional quantum Hall states ^^r— topo- 
logical insulators,— ~— topological superconductors , 30 ! 31 
and two-leg ladders realizing the Haldane gap phase,— ~— 
a real space cut of the system into "halves" will show a 
momentum dependence of the entanglement eigenvalues 
that mimics the physical boundary excitations. In the 
case of a two-leg ladder of spin-1/2 Heisenberg chains 
with a rung coupling; 26 ! 27 the "edge" would be one leg of 
the two. Indeed, tracing out one leg (similar to our trace 
over orbital degrees of freedom in the rung cut) results 
in an entanglement spectrum that closely resembles that 
of the spin-1/2 Heisenbcrg model, which is gaplcssi 26 i 27 

It is interesting to point out that atX = y = 3/4 (see 
FigEJc)), a gapped phase-the non-Haldane phase with 
finite string order ; 49 ' 50 one does not see a characteristic 
spectra of the "edge", which is a spin-1/2 Heisenberg 
chain. This difference compared to the rung coupling 
is intimately related to the very different nature of the 
"four-spin" interaction term (Si ■ Si+x)(fi ■ Ti+i) which 
drives rather different physics between chains i 50 ' 52 

In order to study the competition between the rung 
coupling and the "four-spin" term (Si ■ Si + i)(fi ■ f i+ i) on 
the entanglement spectra in the gapped phase at the ex- 
actly solvable point X = Y = 3/4 of Eq.([T]), we consider 
a perturbation of the form^ 



rung 



J ri 



(23) 



In the spin-ladder interpretation of Hst, this term is 
akin to an AFM exchange coupling between rungs of 
the ladder and hence the reason we have chosen to 
refer to it as rung coupling. This term breaks the 
global SU(2)x SU(2) rotation symmetry down to just 
SU(2). The angular momentum symmetry properties of 
the ground state(GS) are thus modified too. It is no 
longer necessarily true that 5 2 ot = r 2 ot = but only that 
(Stot + ^tot) 2 = 0. This then leads to the entanglement 
spectrum having support in more than one S% ot sector 
of Hs and allows the possibility of higher spin angular 
momentum multiplets in the entanglement spectrum. 

First, we consider the effect of H lung on the entangle- 
ment spectrum of the KM point where our results are 
shown in Fig [12] In Fig fT2T a). we observed that even 
with a very small J run g > relative to X and Y, the 
ES exhibits triplets and higher odd multiplets that de- 
scend from infinite entanglement energy, thus rendering 
A£i ow finite. This ES is significantly different from the 
ES obtained by just perturbing X or Y from the KM 
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FIG. 12. (color online) Momentum resolved entanglement 
spectra for L — 12 chains at the KM point with J 1U ng cou- 
pling (a) J rung = 0.001. (b)J rung = 0.01. (c) J rung = 0.1. (d) 
Jrung = 0.5. In (d) the spectrum resembles the energy spec- 
trum of a Heisenberg spin-1/2 chain and the ES of a two-leg 
ladder without the four spin or dimer-dimer interaction term. 
The degeneracies of the levels are denoted by the symbols 
displayed in the legend. 



point, where only singlet or doublets arc observed due 
to symmetries. However, both share a A£i ow gap which 
indicates that they are in the same phase. It should also 
be noted that when J rung = 0.001 as in Fig fTST a). the 
ground state is still expected to be in the same gapped 
topological non-Haldane phase. This is not as surpris- 
ing as it seems since even with a pure spin-1/2 two 
leg Heisenbcrg ladder without the four spin-interaction 
term, changes in the multiplet structure of the excited 
entanglement energies can occur without going through 
a phase transition . 26 ' 27 In Fig. 3 and Fig. 4 of Ref . [26| . en- 
tanglement spectra were computed for the rung singlet 
(I) and rung singlet (II) phases which exhibit different 
lowest level multiplets (triplets and sextets respectively). 
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Yet, the ground state between these "phases" is a gapped 
product states of rung dimers and is not strictly at a 
phase transition. Hence based on this and our results, 
we conclude that comparing the multiplet structure of 
lowest level ES excitations of a rung cut is not a well 
suited method to distinguish between gapped topologi- 
cal quantum phases. 

Returning to Fig[T2J with increasing Jrung the gap 
A£ low diminishes and more complicated degenerate en- 
tanglement levels appear in the spectrum. The largest 



(a) X=Y=0.25,/ run „=0.15 (b) X=Y=0.25,y mn8 =0.2 



rung coupling value that we considered J rl 



0.5 shown 



in Fig|l2jd) somewhat resembles the degeneracy struc- 
ture of the ES of two leg Heisenberg ladder in the Hal- 
dane and rung singlet (I) phased which mimics the phys- 
ical spectrum of the spin- 1/2 Heisenberg mode h 26 i 27 But 
differences appear in the details of the position of the lev- 
els and the overall dispersion. Nevertheless, in the large 
J rung > limit, we expect to recover a ground state sim- 
ilar to the rung singlet (I) phase and we should recover 
the results of Ref . [26| . 

Finally, we consider perturbing around the SU(4) sym- 
metric point and our results are presented in Fig |13l 
Remarkably the ES is surprisingly robust and remains 
roughly unchanged from the SU(4) point ES even when 
•Aung = 0.15 as shown in FigfT5Ta'). Nevertheless with 
increasing J rung dramatic changes start to appear, in- 
dicative of some sort of phase transition or crossover. 
In Fig fl3r b) with J 1U ng = 0.2, there is a proliferation 
of doublets which occur, even at momentum away from 
p = ±7r/2. There, A£ n igh is still somewhat open. In- 
creasing J rung further causes the A£ n igh to collapse and 
eventually the even multiplcts give way to only odd mul- 
tiplcts in the ES. We observed that the final spectrum 
in Fig|13fd) does not vary much with further increases 
in J rung and the low entanglement energy dispersion ap- 
pears to be rather flat. 

These results on perturbations about the exactly solv- 
able KM point and the SU(4) symmetric point are still 
poorly understood. But they do highlight the rich struc- 
ture that is revealed by the entanglement spectrum when 
changes to the ground state occur, especially in the case 
of perturbing the KM ground state. 



IV. BOND OR LEG CUT ENTANGLEMENT 
SPECTRUM 

In order to obtain complementary information on the 
entanglement properties of the one-dimensional Kugel- 
Khomskii Hamiltonian ([1]), we now consider a cut in 
the Hilbert space where the system is cut in "half 
lengthwise^ through two bonds in real space of the pe- 
riodic chain. This partition then defines two contiguous 
blocks of sites each with length L/2, where one of the 
blocks is traced over. In the spin-ladder interpretation, 
this partition cuts the legs of the ladder symmetrically. 

This is the most common partition used for gapped 
systems, and is natural for revealing the real boundary 
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FIG. 13. (color online) Momentum resolved entanglement 
spectra for L = 12 chains at the SU(4) point with J Tmxs cou- 
pling (a) Jrung = 0.15. (b)Jrung = 0.2. (c) J rung = 0.25. (d) 
Jrung = 0.35. The degeneracies of the levels are denoted by 
the symbols displayed in the legend. 



excitations of an open chain in the Haldanc phase, for 
example^ We study system sizes of L = 8 and 12. In a 
symmetry protected topological phase, like the Haldane 
phase, the lowest levels of the entanglement spectrum will 
be two-fold degenerate per virtual edge due to free spin- 
1/2 edge states ,21 This even degeneracy is protected by 
various symmetries, including inversion symmetry i 24 i 75 
Note that with the periodic boundary conditions we use 
there are always two virtual edges. The quantity of in- 
terest is then the parity of the degeneracies per edge as 
it characterizes the localized edge excitations of the en- 
tanglement Hamiltonian -Kent' 28 ' 76 

We first investigate the entanglement spectrum at the 
KM point. Picking a single MPS groundstate out of 
the two degenerate ones that spontaneously breaks the 
translational symmetry^ we used the analytical formal- 
ism developed in Ref.[24[. We verified that the entan- 
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glement spectrum is four-fold degenerate (two per vir- 
tual boundary) in the thermodynamic limit.— One can 
also intuitively sec this from considering the wavefunc- 
tion which is a staggered dimer covering of a two leg 
ladder. Then a bond or leg cut will always cut a sin- 
glet dimer (which could be a spin or orbital dimer) at 
each edge which leads to the two fold degenerate free 
spin-1/2 levels on each edge. Thus the KM point entan- 
glement spectrum only has one set of two-fold degener- 
ate eigenvalues per virtual edge, exactly as the Afflcck- 
Kennedy-Lieb-TasakiZ£ state. However, the arguments 
developed in Ref.[24[ break down for degenerate ground 
states, which is the case for the KM point. Nevertheless, 
this even degeneracy when properly interpreted is still 
symmetry protected and should remain throughout the 
non-Haldane phase. As was emphasized earlier, gener- 
ically in a finite size system the dimerization is never 
exact and one obtains almost degenerate ground states 
with definite linear momentum K = 0, 7r. In this case 
we find that in the neighborhood of the KM point, there 
are two sets of nearly four-fold degenerate entanglement 
levels as opposed to one. These numerical results are 
presented in Fig ll4l But this may be accounted for by 
taking linear combinations of the MPS states like was 
done in Section IIII B 41 to arrive at definite K ground 
states. Then one can see through linearly superposing 
the individual Schmidt decompositions of the exact MPS 
states that the final reduced density matrix now should 
have two sets of four-fold degenerate entanglement levels 
which is still even per virtual edge. By adiabatic continu- 
ity, if one limits to fixed K ground states these degenera- 
cies should remain invariant as long as the real excitation 
gap remains finite. We emphasize that this added com- 
plication in the entanglement spectrum at finite system 
size is really another artifact of the dimerization order 
exhibited by the non-Haldane phase and a proper inter- 
pretation requires taking into account this degeneracy by 
choosing a definite linear momentum K. 

This special degeneracy can be used to identify the 
topological order in this state, but since it is the similar 
to the Haldane phase the bond cut alone does not al- 
low one to distinguish the Haldane phase from the non- 
Haldane phase. The failure of the degeneracy of the 
entanglement spectrum to identify different topological 
phases has been pointed out before^ (It is possible that 
a "particle partition"— entanglement spectrum might re- 
veal the difference between the Haldane and non-Haldane 
phase as their excitations are different^) However, our 
experience with the rung cut leads us to believe that a 
combination of a leg and a rung cut could provide a way 
to distinguish many different topological gapped phases 
of two leg ladder systems or spin chains with orbital de- 
generacy 

We also study the bond cut entanglement spectrum at 
the SU(4) point. At this point the continuum limit of 
the entanglement spectrum is completely determined' 9 
by the central charge, c = 3. For a finite gapless system, 
the mean number of eigenvalues of the reduced density 
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FIG. 14. (color online) The entanglement spectrum at the 
KM point plotted versus system size. The black dots mark 
the degeneracy of each level. For 8 sites, there is a residual ef- 
fective coupling between virtual edge excitations in H en t that 
splits the four fold degenerate states into singlets and triplets. 
For 12 sites where the edges are more spatially separated the 
4-fold degenerate or 2-fold degenerate per edge is restored to 
a good approximation. The eigenvalues are equal up to the 
third decimal place. 
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FIG. 15. (color online) The mean number of eigenvalues, n(A), 
of the reduced density matrix large r than a given eigenvalue, 
A, plotted as a function of z = ln(A ma a:/A)) in Ea. (|24|) 
for / = 4 and L = 12. The blue circles represent data for 
X = 0.25. The lowest level at z = is non-degenerate and the 
next level is 15-fold degenerate. The green squares represent 
data for X = 0.2. The lowest level at z — is non-degenerate 
and the next level is 9-fold degenerate. The red line is the 
analytical prediction.— 



matrix larger than a given eigenvalue, A , is given by 
n(A) = I (2y/b\ii(X 

/(Mir /A)), (24) 

where Iq is the modified Bessel function, b = 
| ln(-§ sin(7r -|-)), c is the central charge, and / is the sub- 
chain length. For simplicity we let z = 1\Jb ln(A ma;r /A)). 
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In Fig ll5l wc investigate how the low energy excitations 
of (p} effects the degeneracy of the entanglement spec- 
trum in the gapless region V in FigfT] We first compute 
the entanglement spectrum at the SU(4) point. This is 
done a chain of length 12 by tracing over 8 sites. A 
subchain of length 4 is chosen to better reflect the ther- 
modynamic limit. After the lowest entanglement level 
we observe an entanglement level with a 15-fold degen- 
eracy. If we break the SU(4) symmetry by moving along 
the X = Y line to X = 0.2, the symmetry is reduced to 
SU(2) x SU(2). We observe the 15-fold degeneracy splits 
into 9-fold and 6-fold degenerate sets of levels. The en- 
tanglement levels that are 9-fold degenerate have a lower 
entanglement energy and are thus more entangled. We 
can relate these degeneracies of the entanglement spec- 
trum to the degeneracies of the lowest energy excitations. 
At the SU(4) point, excitations with quantum numbers 
(S 2 ,t 2 ) = (0,1), (1,0), (1,1) are degenerate^ (where "1" 
means triplet) and transform under the fundamental ir- 
reducible representation of SU(4) and hence are 15- fold 
degenerate. At X = Y = 0.2, however the lowest energy 
excitation is the (S 2 ,t 2 ) = (1, 1) mode which transform 
as two independent triplets under SU(2)x SU(2). This 
then gives the 9-fold degeneracy. At that point the ex- 
citations with quantum numbers (S 2 ,t 2 ) = (1,0), (0,1) 
are the next lowest excitations and they together are de- 
generate and give the 6-fold degenerate excitations. This 
is somewhat surprising since the ground state is in the 
gapless phase. Unfortunately, our system sizes are too 
small for our data to fall nicely on the curve (|24|) derived 
for the continuum limit. However, both the n(A) and the 
dependence of the entanglement entropy (not shown) on 
sub-system length are consistent with c = 3 given the 
size limitations. 



V. MOMENTUM CUT 

In our final section of the paper, we consider a cut 
defined with momentum space "orbitals" of the Hilbcrt 
space (which is very non-local in position space) at the 
SU(4) point. This momentum-space cut has earlier been 
used to reveal interesting structure in the entanglement 
spectrum, including a clear entanglement gap and CFT 
counting of levels in the spin- 1/2 Hciscnberg chain system 
that reflect the gapless bulk excitations^ In this section, 
we extend this momentum space partition of Ref. [25[ to a 
spin chain with orbital degeneracy or equivalently a two 
leg spin ladder. 

To define this cut, we first write the ground state wave- 
function in a Holstcin-Primakoff representation, 

l*> = 

Y ^(z ji: ...^jjw^, ...,w i JSy i ...SrTr... T r\F), 

ji<—<j K 

ii<...<i K 

(25) 

where \F) = | f • ■ • t) is the fully polarized state and 



K = L/2 is number of spin (S) and pseudospin (r) 
flip operators. K is set by the polarization condition 
5f ot = r t 2 ot = that is satisfied by the ground state. 
We have also expressed the real space wavefunction co- 
efficients in U(l) coordinates of the periodic chain with 
Zj := c l ~^ and Wj := e 1- ^ 2 being positions of a low- 
ered spin and a pseudospin respectively. Note that the 
complex valued function ^{zj t , ...\wi 1 , ...) is totally sym- 
metric in all its arguments and is zero whenever any two 
arguments are equal. This property derives from the mu- 
tual commutativity of these spin flip operators and that 
their square is always zero, {Sj) 2 = {rf ) 2 = 0. Next, 
we define the Fourier and inverse Fourier transforms of 
these spin flip operators by 

L L 
q- — \ " ~™ c- q- _ \ " 7 m 5- 

■ ^jYj ^— ^ 3 3 ' 3 ^/Jj 
l L l L 

^ L 1=1 ^ L rh=l 

(26) 

The ground state wavefunction may then be expressed in 
terms of these Fourier transformed operators as 

l*> = 

^ *(mi, ...,m K |mi, m K )5~ 1 ...S~ \F), 

rhi rh K 

(27) 

where the sum over the different "momenta" rrij and rhi 
takes values in {0, 1, . . . , L — 1} but is otherwise unre- 
stricted. In the rest of this section we will measure the 
crystal momentum in units of 2~k/L. The coefficients of 
([25]) and (f27|) are related to each other according to 

*(mi, ...,m K \rhi, ...,fh K ) = 

Jl<...<7« 
i\ <...<i K 

(28) 

Note that ^ remains symmetric under the action of the 
permutation group of k objects, S K on (mi,...,m re ) and 
(mi, rh K ). However, the spin flip operators and 
in Eq. (|26p may be freely permuted because they are all 
mutually commuting and no longer square to zero unlike 
their real space counterparts. Thus in the sum over m's 
and m's, we should only consider terms that are unique. 
This leads to a description of states which resembles an 
occupation number basis for bosonic particles. Thus in- 
spired by the usual bosonic Fock space, we define a set 
of basis states by 

( C-\n m (~-\n m 

|n,n):= JJ 1 m \ m \ \F), (29) 

m =0 V n m- n m- 
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where n := (no, . . . , Ul-x) and n := (no, . . . , Ul-i)- 
The spin flip operators of ([27]) may be thought of as 
"magnon"-like creation operators acting on the "vac- 
uum" \F) and (n, n) corresponds to an occupation config- 
uration of such magnon orbitals^ However [5 1 ", oc 
and similarly for fjj. Hence, these are not strictly 
speaking bosonic creation operators. Moreover, the basis 
defined by such states is non-orthogonal and overcom- 
plete, although they arc spin, pseudospin and momentum 
cigenstatcs. Nevertheless, we shall refer to the orbital 
states created by S 1 " and as magnons and regard the 
occupation basis as orthonormal. With this in mind, the 
ground state wavefunction may be uniquely written in 
the occupation basis of magnons as follows 

|tf)=^*(n,n)|n,n), (30) 

n.n 

with 

L-l 1 

*(n,n) = TT - y(mi,...,m K \mx,...,m K ), 



where mi < ni2 . . . < m K is a magnon configuration rep- 
resentative of the occupation {no, n£_i} and similarly 
for fh-i < ... < rh K and {fi m , n^-i}. However, these oc- 
cupation numbers are not entirely unrestricted as symme- 
tries of the wavefunction will impose constraints on their 
space of possibilities. For example, the S* ot = r t z ot = 
condition constrains their sum such that 

L-l L-l 

n m = ^2 n m = k. (32) 

m— m— 

Due to time-reversal symmetry, the ground state may 
only have momentum K = or K = L/2 = n in units 
of 2-k/L. This then imposes a constraint on the total 
momentum 

L-l 

M t ot := M A + M B = ^ m x (n m + n m ), 
such that 

M tot mod L = K, (33) 
where K is measured in units of 2it/L. Here, M A '■= 

Em^" 1 ™^"' + and Mb ■= Z)m=L/2 TO ( n i + 

are the total right and left moving momenta from both 
spin(S') and orbital(r) magnons. Finally, we introduce 
the bipartite partition in momentum orbital space by 
splitting the orbitals into right movers with m, m = 
0, L/2 — 1 (partition A) and left movers with m,rh = 
L/2, L — 1 (partition B). The reduced density is then 
formed by tracing out the orbitals in B in the pure ground 
state density matrix, that is /5 re d = TtbI^X^I- 

Next, we discuss the various ways in which we or- 
ganized the the spectrum of H cnt = — In p re( j which is 



the entanglement spectrum. Let Na '■= X)m=o Um 
be the total number of spin magnons in A and Na '■= 
2m=o 1 "™ t ne total number of pseudospin magnons in 
A and similarly for the B partition. Then the sum rule 
of Ea 152]) implies N A + N B = N A + N B = k, which is sat- 
isfied by each ket in the occupation basis expansion. Be- 
cause of this sum rule on the total magnon numbers, p rc( j 
will be block diagonal in (Na,Na) blocks. In our study 
we focus on blocks where the A partition contains half the 
total number of magnons, that is Na = Na = n/2. Thus 
partition A will posses a spin and pseudospin polariza- 
tion of (AS&t, Ar t z ot ) = (-N A , -N a ) relative to \F). In 
addition, we block or resolve the spectrum of Pred|(jv A jv a ) 
into Ma quantum numbers. It must be mentioned, how- 
ever, that Ma is not a strictly a good quantum number as 
was in the case of the spin- 1/2 Heisenbcrg chain.— Only 
when the wavefunction weights (^(rio, -OP have support 
in one M to t sector, can Ma be a good quantum number 
because M tot too is a good quantum number of 5*. Our 
numerical results show a distribution of weights in more 
than one Mtot sector. Nevertheless, in regards to the 
lowest set of levels in the entanglement spectrum, Ma 
remains an approximately good quantum number to la- 
bel these levels. In Appendix [Bj we present additional 
numerical results to substantiate this claim. Hence, each 
Ma sector contains a net spin and pseudospin current of 
(Jf , J*) = (-N A M A , -N A , M a ) in units of 2-k/L. 

We present our numerical computations of the entan- 
glement spectrum at the SU(4) point. Numerical limi- 
tations allow us to only consider systems of sizes L = 8 
and L = 12. Shown in Fig|16l is the momentum cut en- 
tanglement spectrum for a system of size L = 12 at the 
SU(4) point in the sector (N A , N A ) = (3,3). In contrast 
to the a single spin- 1/2 Hcisenberg chain, an entangle- 
ment gap is not exhibited in our spectrum. This goes 
against the expectation that the presence of gapless bulk 
excitations is manifested by an entanglement gap which 
separates a set of universal low lying entanglement ener- 
gies with CFT counting of levels, as was demonstrated 
in the spin-1/2 Heisenberg chain.— Hence, this particular 
partition, as we have performed it, is ill suited to infer 
the existence of gapless bulk excitations from the entan- 
glement spectrum in this model. The complications with 
the several M to t sectors may have been a precursor of 
this fact. Nevertheless, we next consider the counting 
of low energy entanglement levels when they are well- 
defined. Due to the periodic nature of momentum space 
(first Brillouin zone) , there are really two cut regions be- 
tween m = L — 1,0 and between m = L/2 — l,L/2. 
Thus we should be expect to see CFT counting of levels 
from either low or high momentum. However, we observe 
instead a counting of lowest energy levels which is asym- 
metric from high and low Ma, which again is in contrast 
with the single Heisenberg chain. From large Ma we see 
a count of lowest energy levels of 1,2,5,10. From small 
Ma where there appears to be a small gap opening we 
observe a count of 1,2,5,8 instead, but this may be due to 
finite size effects as it was not observed previously with 
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FIG. 16. The momentum cut entanglement spectrum at the SU(4) symmetric point for L = 12 in the sector (Na,Na) = (3, 3). 
The ground state possesses momentum K = k = L/2 in units of 2n/L. The count of the lowest levels from high Ma (right) is 
1,2,5,10 and from low Ma (left) is 1,2,5,8. The floating numbers in the plot denote the total number of levels directly below. 
There are two pairs of levels which are almost degenerate at Ma = 27 which are not resolved at this scale. The spectrum does 
not exhibit an entanglement gap despite being in the gapless phase V described by a sum of two SU(2)2 WZW theories or a 
single SU(4)i WZW theory. 



L = 8 data (See Fig flU)) . Moreover, the lowest energy 
levels of the spectrum also appears to exhibit a strong 
dispersion which is be bowl-like in contrast to the rela- 
tively fiat set of lowest levels seen in the spin- 1/2 Heisen- 
berg chain. We note that the 1,2,5,10 sequence may be 
obtained from the dimensions of highest weight and de- 
scendant states of a tensor product of two independent 
Verma modules of the Virasoro algebra, each with the 
usual counting 1,1,2,3 with increasing level. Presently, 
we lack satisfactory interpretations of the absence of the 
entanglement gap and the observed counting of states 
in terms of the CFT that describes the low energy bulk 
excitations. 



VI. CONCLUSIONS 

In conclusion, we have investigated the entanglement 
properties of the one-dimensional Kugel-Khomskii model 
([T]) using exact diagonalization for system sizes of L = 
8, 12 and 16 sites. By considering various partitions or 
"cuts" of the Hilbert space, including a rung (orbital), 
bond, and momentum cut we have produced a compre- 
hensive picture of the phase diagram from the point-of- 
view of the quantum entanglement. 

Of the three partitions that we have studied, the par- 
tition into spin and orbital degrees of freedom was the 
most intensely studied. With that partition, we have 
shown that there are crucial finite-size effects at L = 8 
in the entanglement entropy and the entanglement spec- 
trum that impacted earlier exact diagonalization stud- 
ies of this model which led to a misinterpretation of the 



phase diagram^ In particular, using the entanglement, 
we are able to reaffirm that regions IV and VI are adi- 
abatically connected to one another once ground state 
degeneracies due to dimerization are taken into account. 
We have observed and explained various symmetries of 
the entanglement spectrum obtained by this cut which 
can depend crucially on the ground state linear momen- 
tum K . More importantly we observed that the entan- 
glement spectrum of the rung cut does not in any way 
resemble the energy spectrum of a single spin-1/2 Heisen- 
berg chain which is the system one obtains when the or- 
bital and spin degrees of freedom are decoupled. Thus 
the Kugel-Khomskii chain offers a counter example to 
the commonly held boundary-bulk correspondence con- 
jecture of the entanglement spectrum*^— This is inti- 
mately related to the natural c = 3 CFT of the gap- 
less phases, while two spin-1/2 Heisenberg edges only 
produce a net central charge of c = 2. With this cut, 
we have also observed two important entanglement gaps. 
The lower gap which we have denoted as A£i ow charac- 
terizes the dimerization of phase IV, which is the non- 
Haldane phase. A£i ow separates the two lowest entangle- 
ment levels from a continuum and diverges at the exactly 
solvable Kolezhuk-Mikeska (KM) point. We have com- 
puted the entanglement spectrum exactly at this special 
point. The second and higher gap which we have de- 
noted by A£high characterizes the gapless AFM phase or 
phase V. It separates a set of N(L) = 2 L / 2 ~ 1 levels from 
the continuum and we show through studies of the cor- 
relation functions that these levels are the source of the 
critical fluctuations of the ground state spin-spin corre- 
lators. Our study of the evolution of the entanglement 
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spectra through the Kosterlitz-Thouless phase transition 
between phase V and IV (VI) sharply identifies the level 
crossings between the approximately degenerate ground 
states. More importantly, once K is fixed as was done in 
Fig Etc), the evolution of the spectrum reveals the open- 
ing and closures of the entanglement gaps which is indica- 
tive of a phase transition. However, our finite size data 
is still unable to definitively identify the SU(4) point as 
the exact location of the critical point of the KT tran- 
sition. By perturbing Eq.([T|) with a direct spin-orbital 
exchange perturbation, we showed that the SU(4) point 
rung cut entanglement spectrum and the A£high gap is 
surprisingly robust. By contrast, perturbations at the 
KM point show a finite A£i ow gap developing but reveal 
odd multiplets in the entanglement levels above A£i ow 
even when remaining in the non-Haldane phase. Thus 
we showed that the multiplct structure of the excited en- 
tanglement levels are a poor identifier of a topological 
phase with this cut. 

With the bond or leg cut entanglement spectrum 
study we observed in the non-Haldane phase two 2-fold 
degenerate entanglement levels per virtual edge which 
we attribute to nearly degenerate ground states due to 
spontaneous dimerization. In the gaplcss phase, we see a 
distribution of eigenvalues which because of small system 
sizes does not fit accurately to CFT predictions^ but is 
still consistent. Furthermore, the lowest entanglement 
level degeneracies reveal the degeneracies of lowest real 
energy excitations. 

Finally, we extended the non-local momentum cut of 
Thomale et. al<^ to two leg ladders or spin chains with 
orbital degeneracy. We obtained the entanglement spec- 
trum at the SU(4) symmetry point for a chain of size 
L = 12. In contrast to the single spin-1/2 Hcisenberg 
chain we do not see a clear entanglement gap which would 
be an indicator of gapless CFT bulk excitations. How- 
ever we observed a new counting of levels that has yet to 
be identified with the SU(4)i or two SU(2) 2 WZW theo- 
ries that are the two alternate low energy descriptions of 
the system. 
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Appendix A: Features of the Entanglement 
Spectrum when K = n 

Our data on the entanglement spectrum reveals non- 
trivial relations between spectra and possible degenera- 
cies when the ground state has momentum K = tt. In 
this appendix we present arguments to explain the origin 
of these features. 

First consider a non-degenerate ground state with 
generic X ^ Y with K = 0, tt and Schmidt decompo- 
sition, 

P <Pp 

and energy H(X,Y)\^>) = E \ty). Then because 
H<y,X)F\$) = TH{X,Y)\$) = E F\V), thus F\<S>) is 
a ground state of H(Y, X) which is also non-degenerate. 
But 

p <f> P 

= EE c-^-- P ) /2 | P 0^)®i(if-p)0 p ), 

P <P(K-p) 

(A2) 

where in the last line we have relabeled the dummy 
momentum p — > K — p. Thus J-\*5f) yields an entan- 
glement spectrum of {(p, £^ (Jf _ p) )}. This implies that 
the (non-degenerate) ground states related by T have 
their entanglement spectra shifted according to K — p 
relative to each other. Since K = or tt, this means 
that |\&) with spectrum {(p, ^0 p )} is related to the spec- 
trum of F\ "J) with spectrum {(p, £<p_ p )} for K = or 
{(p, 6^-p)} for K = tt. Since we have already estab- 
lished the {(p,£<p p )} = {(p, symmetry due to in- 
version I, the K = case requires that the spectra of 
and |^) be identical. But for K = tt the shift is 
non-trivial and is confirmed by our numerics as shown in 
FigEta-b). 

Next focusing on K = tt and the special momenta 
where tt — p = p mod 2tt. This occurs when p = ±7r/2 
and at these special momenta, the |^) and J - ^) levels 
are identical. 

I*} 

{(■J J £0^/2)} = {(2' £0*72)} f A q\ 

{(-|)^_„/ 2 )} = {( - f>^_,/ 2 )} 

Finally, we consider the special limit when X = Y. 
In this case H{X,X) does commute with T and non- 
degenerate ground states transform into themselves mod- 
ulo a sign under T . Considering still only the case when 
K = tt and so we have an additional symmetry in the 
entanglement spectrum, 

{(P, = {(t - P, &,)} = {(■* + P, &_,)} 

= {(*+P,S* f )}- (A4) 
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Hence, p and p + it have identical spectra too and in 
particular p — and p = n are identical. When p = ±7r/2 
we have the empty statement that 



{(±7r/2,^ ±7r/2 )} = {(±7r/2,6^ /2 )} 



(A5) 



For definiteness, take p = it/2 (an identical consider ap- 
plies for p = —it/ 2) and looking at the contribution to 
\ty) in the Schmidt decomposition, 



n \n/2 fa /3 ) <g> \n/2 ff/2 ) 



(A6) 



07T/2 

Applying J 7 yields 

e" C *" /2 K/2 ® |tt/2 V2 ) (A7) 

<K/2 

which is nothing but the statement that spectrum at 
p = it/2 maps back into itself. However if we suppose 
J 7 !*) = — then we cannot have (j>„/ 2 K e l se we 
will arrive at a contradiction (this however is consistent 
with J 7 !*) = 1^))- Hence we must have that 4 > tt/2-L4 i w/2 
and that the eigenvalue , 2 must be degenerate with 
£^ — -. Thus, the entanglement energies must be at least 

two- fold degenerate at p = tt/2 in this case. This then 
explains the doublets seen in Fig|6jd). 



Appendix B: Ma as an approximately good quantum 
number in the momentum cut 



In this appendix we show that blocking p rc d\^ A jy A ) 
into total right moving momentum Ma sectors has little 
quantitative influence on lowest levels of the full entan- 
glement spectrum of p r ed\(M A n a )- I n the case of a single 
spin-1/2 Heiscnberg chain, the expression of ground state 
wavefunction in terms of the magnon occupation basis 
yielded coefficients with weig hts |*(n ,...)| 2 that were 
sharply concentrated around sector Mtot = k 2 = L 2 /4.— 
This led to a sum rule Ma + Mb = k which in turn led to 
/9 ro d being block diagonal in Ma- Unfortunately, in the 
case of the spin-orbital chain, things are more compli- 
cated and the weights are distributed over three to four 
sectors depending on K . When K = the weights are 
distributed over M tot = 2k 2 , 2k 2 ± L and when K = L/2 
they are distributed over M tot = 2k 2 ± k, 2k 2 ± k ± L. 
Shown in FigfT7| is the weight distribution of the coeffi- 
cients for an L = 8 Kugel-Khomskii spin chain ground 
state at the SU(4) symmetric point X = Y = 1/4 where 
most of the weight is concentrated at 2k 2 . 

Nevertheless, it appears that M to t is an almost a good 
quantum number for |\&), at least with regards to charac- 
terizing the entanglement spectrum of p ro d = Ti's|5')( , I'|. 
If we compare the spectrum obtained by just diagonal- 
izing p rc d and the spectrum obtained by first explicitly 



blocking into Ma sectors (ignoring off diagonal terms be- 
tween sectors) and then computing the resulting spec- 
trum, we find spectra which agree very well with one 
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FIG. 17. The total weight of the occupation basis coefficients 
taken from the ground state at the SU(4) point. The system 
size is L — 8 and the ground state momentum is K = 0. 
The distribution of weights in distributed across three sectors, 
Mtot = 2k 2 , 2k 2 ± L, where k = L/2 = 4. 
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FIG. 18. Entanglement spectra arranged in ascending order 
taken from the SU(4) symmetric point of an L — 8 chain 
in the sector (Na,Na) = (2,2). Shown by the boxes (□) is 
the spectrum from first blocking p rc d and the crosses (+) the 
spectrum without blocking p rc d- The two spectra agree at the 
lowest energies and differ at the highest levels £ > 40. 



another. Shown in Fig |18l are the two spectra computed 
with these two different methods. The data shows an 
agreement between the entanglement energies up till the 
very highest levels around £ > 40, where we believe the 
spectrum is less likely be of strong physical significance. 
We also present in Figfl8] the entanglement spectrum re- 
solved in Ma sectors for an L = 8 chain. 

With a system of size L = 12, the agreement between 
blocked and unblocked spectra is not as strong but the 
lowest entanglement energy levels nevertheless are inde- 
pendent of the Ma blocking with significant deviations 
only starting to appear around £ ~ 12. Moreover, we 
expect that only these robust lowest entanglement levels 
will reflect the low real energy properties of the system. 
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FIG. 19. The momentum cut entanglement spectrum of the 
Kugel-Khomskii chain at the SU(4) symmetric point for L — 
8. The spectrum is taken from into Ma blocks of p rc d in the 
sector (Na, Na) = (2, 2). The count of the lowest levels from 
high Ma (right) is 1,2,5 and from low Ma (left) it is 1,1,3. 
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